Giant electro-thermal conductivity and spin-phonon coupling in an antiferromagnetic 

oxide 
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The application of weak electric fields (< 100 V/cm) is found to dramatically enhance the lattice 
thermal conductivity of the antiferromagnetic (AF) insulator CaMnOa over a broad range of temper- 
ature about the Neel ordering point (125 K). The effect is coincident with field-induced de-trapping 
of bound electrons, suggesting that phonon scattering associated with short- and long-ranged AF 
order is suppressed in the presence of the mobilized charge. This interplay between bound charge 
and spin-phonon coupling might allow for the reversible control of spin fluctuations using weak 
external fields. 

PACS numbers: 66.70.-f,75.47.Lx, 75.80.-l-q, 71.55.-i 



The strong magneto-elastic coupling present in cer- 
tain antiferromagnetic (AF) ferroelectrics [l|] (e.g., hexag- 
onal YMnOa) has attracted considerable attention re- 
cently because of its importance in mediating a cou- 
pUng between magnetic and electric orders 0] and the 
potential of such multiferroics for applications 0]. The 
hexagonal manganites also possess very strong spin fluc- 
tuations characteristic of a spin liquid well above their 
Neel temperatures, possibly related to geometric spin 
frustration. This raises the prospect of manipulating 
dynamic magneto-elastic coupling in the paramagnetic 
phase (PM) of these or other compounds using external 
fields. 

Dramatic signatures of such coupling in YMnOs are 
a suppressed thermal conductivity 0] (k) and anoma- 
lous ultrasonic attenuation over a broad tempera- 
ture range in the PM phase. Two AF oxides with quite 
similar features in their thermal conductivities are the 
transition-metal monoxide, MnO Q (T^r = 118 K), 
and orthorhombic CaMnOs 0, [| {Tn ~ 125 K). Like 
YMnOa, both of these materials possess large exchange 
striction [13] (with changes in lattice constants at 
T < Tn, Aa/a > lO"'^), significant spin frustration 
[H [H e/TAT - 4 - 5 (6 is the Curie- Weiss tempera- 
ture) , and strong short-range spin correlations extending 



well above Tn [III, ll3| . The suppressed thermal conduc- 
tivities in the PM phase of YMnOa and CaMnOa (and 
by extension, MnO) have been proposed to arise from 
the scattering of acoustic phonons by nanoscale strains 
generated by short-ranged spin correlations [3,[7|. This 
scattering diminishes rapidly when long range order is es- 
tablished, giving rise to sharp increases in k at T < T^- 
CaMnOs is distinguished from YMnOs and MnO by 
the presence of low-lying electron donor levels (oxygen 
vacancies) from which a small density of mobile electrons 
(An - 10 



16 (^j^ 3 -J ^g^j^ |-,g released under relatively weak 



applied electric fields (F < 100 V/cm) [ij]. Here we 
demonstrate that these mobilized carriers are associated 
with a substantial reduction in phonon scattering, yield- 
ing relative changes in k with field, (1/k)(Ak/AF) that 



are two orders of magnitude larger than found in quan- 
tum paraelectrics such as KTaOa and SrTiOs jlSl where 
the field couples to soft-mode phonons. The apparent ab- 
sence of permanent electric dipoles or lattice anomalies 
in CaMnOa leads to the hypothesis that the mobilized 
electrons themselves mediate a suppression of phonon 
scattering from strain induced by short- and long-ranged 
magnetic order. 

Measurements were performed on single- crystal and 
polycrystalline specimens of CaMnOa. Their prepara- 
tion methods and physical properties have been reported 
elsewhere [l,[li, 16, 17, 3|. We focus here on the crystal 
(dimensions 2 x 0.9 x 0.05 mm'^) for which k was mea- 
sured for two different oxygen configurations produced 
by annealing at 600° C in flowing oxygen and air, respec- 
tively. The net carrier density for the oxygen-annealed 
state, as determined from room-temperature Hall mea- 
surements 55], was ~ 6 X 10^^ cm~'^. This corre- 
sponds to '--^ 3 X 10~* electrons per formula unit, at- 
tributed to a very small oxygen deficiency. The carrier 
density after annealing in air is estimated [3] from the 
increase in room-temperature thermopower (-530 /^/K 
to -370 n/K) to be Af ~ 4 X 10^^ cm"^. k was mea- 
sured with a steady-state technique employing a heater 
and 25 //m-diameter chromel-constantan thermocouple; 
four-probe electrical resistivity (p) was measured during 
the same experiments. k(J, T) was measured for the air- 
annealed crystal in the presence of dc transport currents 
(/ = 5 nA-10 mA). With the specimen suspended in 
vacuum and thermally anchored at only one end. Joule 
heating was significant at the highest currents. The aver- 
age specimen temperature (relative to a resistance sensor 
mounted on the cold stage) was monitored with a second 
thermocouple following application of a transport cur- 
rent, and allowed to stabilize prior to energizing the spec- 
imen heater for determining k. Linearity in the heater 
power - AT response was confirmed at various tempera- 
tures throughout the measurement range. 

Figure[T] (lower panel) shows k{T) after oxygen (dashed 
curve) and air (circles) annealing. The / = and 
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FIG. 1: (Color online) Lower panel (left ordinate): n{T) for 
the air-annealed CaMnOa crystal (circles) at various trans- 
port currents, from right to left (in mA): 10, 5, 3, 2, 1, 0.5, 
0.25, 0.1, 0.05, 0.01, 0. The solid curves bridging k,{T) plots 
are contours of constant electric field (values labeled). The 
dashed curve represents data for the oxygen- annealed crystal 
for / = 0. Upper panel (right ordinate): p{T) for the same 
currents as ii{T) as well as the following (in ^A): 5, 3, 1, 0.3, 
0.1, 0.03, 0.01, 0.005. Inset: current dependence of the slope, 
dK/dT, for air annealing. 



/ = 10 niA plots for air annealing are labeled. Con- 
stant electric field (F) contours (with values labeled) are 
represented by solid curves bridging the constant current 
plots. These fields, F = pj ^ were determined from p{T) 
data (upper panel). The abrupt upturn in k upon enter- 
ing the AF ordered phase is evident near Tjv ~ 120 K. 
Note that for the J = curves the slope just below T/v, 
(i/«/dT|Tjv, is greater for oxygen annealing, and increases 
with / for the air annealed crystal (inset. Fig. [T]). These 
features are discussed further below. The data were re- 
producible upon thermal and current cycling, and after 
subsequent re-annealing. 

Figured (a) shows the T-dependence of the fractional 
change in thermal conductivity, (Ak/k) = [k{F,T) — 
k(0, r)]/K(0, T), computed for these values of F. The en- 
hancement of K, in field is substantial, approaching 100% 
for fields F > 100 V/cm in the range 80 K < T < 130 K, 
and becomes negligible at lower temperatures. Quali- 
tatively similar results were found for the polycrystalline 
specimen, with a temperature interval for the largest field 
effects lower by about 50 K. 

The field-enhanced k correlates with an increased elec- 
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FIG. 2: (Color online) (a) Fractional change in thermal con- 
ductivity vs temperature at constant electric fields, (b) Field- 
induced change in electrical conductivity for the same field 
values as in (a). 

trical conductivity. Act = a{F, T) - cr(0, T) [Fig. [2] (b)]. 
As demonstrated by prior Hall effect studies on this same 
crystal [l3| , the mechanism underlying the increase in a is 
Poole-Frenkel [l^l release of carriers from Coulomb traps 
(oxygen vacancy donor levels) via field-assisted barrier- 
lowering. However, it is clear that the thermal conduc- 
tivity of these mobilized electrons cannot account for 
the increase in k. For example. Act ~ 0.5 Q^^mr~^ at 
F = 144 V/cm and T = 126 K, corresponding to an 
increase in the electron density. An ~ 3 x 10^^ cm~^ 
(using a mobility [11] ^ 1 cm^/V-s), about 50% of the 
equilibrium value at that temperature. The Wiedemann- 
Franz law, which provides an upper-bound estimate 
on the electronic thermal conductivity, yields An,, — 
LqActT ~ 1.5 X 10^^ W/niK, nearly seven orders of 
magnitude smaller than observed. Spin-wave heat con- 
duction, though a possible contributor to k at low T in 
the ordered phase, cannot account for k well above T/v. 
Thus the increase of k must have its origin in reduced 
phonon scattering. 

Bound electrons may cause phonon scattering through 
static local distortions of bond lengths 2l| or more ex- 
tended, dynamic polaronic distortions, and possibly the 
diminution of this scattering as electrons become mobile 
in field explains the observations. But this proposal must 
be viewed as implausible given that the fraction of bound 



electrons released is very small, ^ An/N < 10~^, so 
that the density of such distortions is altered negligibly in 
field. Further supporting this view, k at the highest F for 
the air-annealed crystal exceeds that in zero-field for the 
oxygen-annealed state (Fig. [T|), in spite of the fact that 
the bound electron density in the former exceeds that of 
the latter by an order of magnitude. In addition, there 
is no evidence of permanent electric dipoles in this mate- 
rial 
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122] , the alignment of which in weak fields might 
conceivably alter the phonon scattering. This leads us to 
consider that the mobilized electrons themselves mediate 
a suppression of phonon scattering. 

In order to assess how this might occur and to gain 
insight into the origin of the strong PM-phase phonon 
scattering, it is instructive to examine changes in struc- 
ture induced by the onset of magnetic order. Anoma- 
lous Mn motion, stabilized at T < Tjv, was recently pro- 
posed to explain the k behavior of hexagonal YMnOa 
[23|. However, the unusual Mn-0 bondlength changes [l| 
motivating this proposal are absent in the average and 
local structure of CaMnOa 2J]. On the other hand, ro- 



tations of the MnOg octahedra couple strongly to the 
magnetic order in CaMnOs 
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and provide a mech- 
anism through which spin fluctuations at T > Tn may 
induce strain that scatters heat-carrying acoustic modes. 
Perturbations due to the polaronic motion 14, [l^ 25 1 
of mobilized electrons may compete with the distortions 
favored by short-ranged AF ordered domains. 

Consider the scenario for phonon scattering by spin 
fluctuations introduced by Sharma et al. Strain 
fields associated with short-ranged AF spin-ordered re- 
gions have spatial extent given by the spin correlation 
length ^. Slowly fluctuating on the timescale of lat- 
tice vibrations, they scatter acoustic phonons. For the 
case Ap <§; ^ (where Xp is the phonon wavelength [1^), 
their scattering rate can be approximated (for a spher- 
ical spin-correlated domain) as, t"^^ = povit{£^/2)'^ (pq 
and V are the density of scatterers and phonon velocity, 
respectively). This scattering should be sharply dimin- 
ished at T < T/v as long-range magnetic order is estab- 
lished. Residual strain in AF domain walls and other 
defects (e.g. twins) likely scatter phonons to low tem- 
peratures as evidenced by a phonon mean- free path [i^] , 
'^ph = Sk/Cv^ that remains < 1/im at the lowest mea- 
sured T. 

To examine this picture more quantitatively we em- 
ploy the Debye- Calloway model [28| to compute the lat- 
tice thermal conductivity, adding r^^g to a sum of tradi- 
tional phonon scattering rates arising from other phonons 
(Umklapp), dislocations, and point-like defects [29ii]. We 
restrict our analysis to the PM phase and for simplic- 
ity ignore possible temperature dependencies of no and 

taking 7 = pon{^/2)'^ as a constant parameter for 
T > Tpf. The 7 = data for air and oxygen anneal- 
ing were first fitted by assuming the traditional phonon 
scattering terms to be independent of the annealing treat- 
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FIG. 3: (Color online) Comparison of calculated (curves) and 
measured (symbols) k{T) for 1 = 0. The dashed curve is the 
computed k in the absence of magnetic scattering (7 = 0). 
The solid curves used 7 = 2.51(0.82) x lO'^m"^ for air (oxygen) 
annealing. Inset: magnetic scattering strength versus field- 
induced carrier density determined from fits to k{F) at fixed 
T. 



ment, and the computed k in the absence of magnetic 
scattering (7 = 0) was constrained to match the oxygen 
annealed data at the lowest T and to exceed the high- 
est measured values (including those for air annealing in 
applied field) at intermediate T. Subsequently, the data 
for air annealing in applied field were fitted with 7 as the 
only adjustable parameter. 

Fig. [3] compares the calculated (solid curves) and mea- 
sured k(T). Also shown is the hypothetical k{T) in the 
absence of magnetic scattering (7 = 0, dashed curve). 
The values of 7 = 2.51(0.82) x lO'^ m'^ for air (oxy- 
gen) annealing are similar to that employed for YMnOa 
[4|. For self-consistency, the distance between scatter- 
ers d = 2(3/47rno)^^'^, must be larger than their size ^. 
Both YMnOs 0] and MnO 11 1 have comparable values 
of ^ ~ 20 — 50A. In the absence of a direct measure- 
ment of ^ for CaMnOs, taking ^ = 40A implies a density 
of scatterers po — 2.0 x 10^^ cm~'^ for air annealing, or 
equivalently, d ~ lOOA > ^ as required. 

Fitting the model to k{F) at fixed T yields the decrease 
in magnetic scattering vs. the relative change in mobile 
carrier density (inset of Fig. The decrease in 7, by 
up to a factor of 10 at the highest applied fields near T/v, 
implies a decrease in poj C or both. 

The nucleation of AF domains in the PM phase and 
their growth in the ordered phase likely depend on 
domain- wall pinning near lattice defects [30[. Mobilized 
electrons could reduce po by depinning AF domains, an 
idea motivated by comparing the zero-field data for air 
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and oxygen annealing. The larger k and dK/dT\Tfj ob- 
served for the more oxygenated crystal in the ordered 
phase (Fig. [T]) suggests a coarsening of the domain struc- 
ture, consistent with observations in Gd2(Mo04)3 of 
stron g p honon scattering at ferroelastic domain bound- 



3l|. Oxygen vacancies or their complexes are thus jgj 



aries 

implicated in pinning (the estimate of po above suggests 
that ~ 10% of vacancies participate in pinning). 

The correlation between k and dK/dT\x^ in applied 
electric field (Fig. [T|) suggests that mobihzed electrons 
cause depinning. The polaronic motion of charge carri- 
ers can increase the frequency of spin fluctuations [s^ . 
effectively reducing pq. Regarding ^ there is prece- 
dent for its diminution by light charge-carrier doping in 
Cai_2:Laa;Mn03 [1^ and La2-a;Sra;Cu04 [32]. However, 
the carrier densities in these cases are several orders of 
magnitude larger than An induced by field in the present 
experiments. Both n and dn/dTlx,^ decrease with sub- 
stitutional doping in Cai_3;Laa,Mn03 Q, contrary to the 
field-induced trend, but this may simply reflect the pre- 
dominance of additional scattering from substitutional 
disorder. More direct measurements of the AF correla- 
tions in CaMnOa as functions of temperature and applied 
electric fleld are clearly called for to further test and re- 
fine the ideas put forward here. 

In summary, the present measurements of thermal con- 
ductivity in CaMnOa imply that a small density of elec- 
trons released from donor states in weak electric fields 
substantially reduce the scattering of phonons associ- 
ated with short- and long-ranged AF order. The re- 
sults suggest that the correlated spin volume in the PM 
phase is highly sensitive to the injected carrier density 
and reversibly controlled by a transport current. Strong 
AF spin fluctuations, magnetostriction, and the presence 
of loosely-bound charges in defect states appear to be 
key ingredients underlying the observations. This novel 
manipulation of dynamic magneto-elastic coupling mo- 
tivates further studies of spin correlations in this com- 
pound or others having similar characteristics. 
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